Like/dislike analysis using EEG: determination of most discriminative channels and frequencies.

In this study, we have analyzed electroencephalography (EEG) signals to investigate the following issues, (i) which frequencies and EEG channels could be relatively better indicators of preference (like or dislike decisions) of consumer products, (ii) timing characteristic of "like" decisions during such mental processes. For this purpose, we have obtained multichannel EEG recordings from 15 subjects, during total of 16 epochs of 10 s long, while they were presented with some shoe photographs. When they liked a specific shoe, they pressed on a button and marked the time of this activity and the particular epoch was labeled as a LIKE case. No button press meant that the subject did not like the particular shoe that was displayed and corresponding epoch designated as a DISLIKE case. After preprocessing, power spectral density (PSD) of EEG data was estimated at different frequencies (4, 5, …, 40 Hz) using the Burg method, for each epoch corresponding to one shoe presentation. Each subject's data consisted of normalized PSD values (NPVs) from all LIKE and DISLIKE cases/epochs coming from all 19 EEG channels. In order to determine the most discriminative frequencies and channels, we have utilized logistic regression, where LIKE/DISLIKE status was used as a categorical (binary) response variable and corresponding NPVs were the continuously valued input variables or predictors. We observed that when all the NPVs (total of 37) are used as predictors, the regression problem was becoming ill-posed due to large number of predictors (compared to the number of samples) and high correlation among predictors. To circumvent this issue, we have divided the frequency band into low frequency (LF) 4-19 Hz and high frequency (HF) 20-40 Hz bands and analyzed the influence of the NPV in these bands separately. Then, using the p-values that indicate how significantly estimated predictor weights are different than zero, we have determined the NPVs and channels that are more influential in determining the outcome, i.e., like/dislike decision. In the LF band, 4 and 5 Hz were found to be the most discriminative frequencies (MDFs). In the HF band, none of the frequencies seemed offer significant information. When both male and female data was used, in the LF band, a frontal channel on the left (F7-A1) and a temporal channel on the right (T6-A2) were found to be the most discriminative channels (MDCs). In the HF band, MDCs were central (Cz-A1) and occipital on the left (O1-A1) channels. The results of like timings suggest that male and female behavior for this set of stimulant images were similar.

Analysis of coronary angiography related psychophysiological responses.

BACKGROUND: Coronary angiography is an important tool in diagnosis of cardiovascular diseases. However, it is the administration is relatively stressful and emotionally traumatic for the subjects. The aim of this study is to evaluate psychophysiological responses induced by the coronary angiography instead of subjective methods such as a questionnaire. We have also evaluated the influence of the tranquilizer on the psychophysiological responses. METHODS: Electrocardiography (ECG), Blood Volume Pulse (BVP), and Galvanic Skin Response (GSR) of 34 patients who underwent coronary angiography operation were recorded. Recordings were done at three phases: "1 hour before," "during," and "1 hour after" the coronary angiography test. Total of 5 features obtained from the physiological signals were compared across these three phases. Sixteen of the patients were administered 5 mg of a tranquilizer (Diazepam) before the operation and remaining 18 were not. RESULTS: Our results indicate that there is a strong correlation between features (LF/HF, Bk, DN1/DN2, skin conductance level and seg_mean) in terms of reflecting psychophysiological responses. However only DN1/DN2 feature has statistically significant differences between angiography phases (for diazepam: p = 0.0201, for non_diazepam p = 0.0224). We also note that there are statistically significant differences between the diazepam and non-diazepam groups for seg_mean features in "before", "during" and "after" phases (p = 0.0156, 0.0282, and 0.0443, respectively). CONCLUSIONS: The most intense sympathetic activity is observed in the "during" angiography phase for both of the groups. The obtained features can be used in some clinical studies where generation of the customized/individual diagnoses styles and quantitative evaluation of psychophysiological responses is necessary.

Classification of root canal microorganisms using electronic-nose and discriminant analysis.

BACKGROUND: Root canal treatment is a debridement process which disrupts and removes entire microorganisms from the root canal system. Identification of microorganisms may help clinicians decide on treatment alternatives such as using different irrigants, intracanal medicaments and antibiotics. However, the difficulty in cultivation and the complexity in isolation of predominant anaerobic microorganisms make clinicians resort to empirical medical treatments. For this reason, identification of microorganisms is not a routinely used procedure in root canal treatment. In this study, we aimed at classifying 7 different standard microorganism strains which are frequently seen in root canal infections, using odor data collected using an electronic nose instrument. METHOD: Our microorganism odor data set consisted of 5 repeated samples from 7 different classes at 4 concentration levels. For each concentration, 35 samples were classified using 3 different discriminant analysis methods. In order to determine an optimal setting for using electronic-nose in such an application, we have tried 3 different approaches in evaluating sensor responses. Moreover, we have used 3 different sensor baseline values in normalizing sensor responses. Since the number of sensors is relatively large compared to sample size, we have also investigated the influence of two different dimension reduction methods on classification performance. RESULTS: We have found that quadratic type discriminant analysis outperforms other varieties of this method. We have also observed that classification performance decreases as the concentration decreases. Among different baseline values used for pre-processing the sensor responses, the model where the minimum values of sensor readings in the sample were accepted as the baseline yields better classification performance. Corresponding to this optimal choice of baseline value, we have noted that among different sensor response model and feature reduction method combinations, the difference model with standard deviation based dimension reduction or normalized fractional difference model with principal component analysis based dimension reduction results in the best overall performance across different concentrations. CONCLUSION: Our results reveal that the electronic nose technology is a promising and convenient alternative for classifying microorganisms that cause root canal infections. With our comprehensive approach, we have also determined optimal settings to obtain higher classification performance using this technology and discriminant analysis.

Sleep stage and obstructive apneaic epoch classification using single-lead ECG.

BACKGROUND: Polysomnography (PSG) is used to define physiological sleep and different physiological sleep stages, to assess sleep quality and diagnose many types of sleep disorders such as obstructive sleep apnea. However, PSG requires not only the connection of various sensors and electrodes to the subject but also spending the night in a bed that is different from the subject's own bed. This study is designed to investigate the feasibility of automatic classification of sleep stages and obstructive apneaic epochs using only the features derived from a single-lead electrocardiography (ECG) signal. METHODS: For this purpose, PSG recordings (ECG included) were obtained during the night's sleep (mean duration 7 hours) of 17 subjects (5 men) with ages between 26 and 67. Based on these recordings, sleep experts performed sleep scoring for each subject. This study consisted of the following steps: (1) Visual inspection of ECG data corresponding to each 30-second epoch, and selection of epochs with relatively clean signals, (2) beat-to-beat interval (RR interval) computation using an R-peak detection algorithm, (3) feature extraction from RR interval values, and (4) classification of sleep stages (or obstructive apneaic periods) using one-versus-rest approach. The features used in the study were the median value, the difference between the 75 and 25 percentile values, and mean absolute deviations of the RR intervals computed for each epoch. The k-nearest-neighbor (kNN), quadratic discriminant analysis (QDA), and support vector machines (SVM) methods were used as the classification tools. In the testing procedure 10-fold cross-validation was employed. RESULTS: QDA and SVM performed similarly well and significantly better than kNN for both sleep stage and apneaic epoch classification studies. The classification accuracy rates were between 80 and 90% for the stages other than non-rapid-eye-movement stage 2. The accuracies were 60 or 70% for that specific stage. In five obstructive sleep apnea (OSA) patients, the accurate apneaic epoch detection rates were over 89% for QDA and SVM. CONCLUSION: This study, in general, showed that RR-interval based classification, which requires only single-lead ECG, is feasible for sleep stage and apneaic epoch determination and can pave the road for a simple automatic classification system suitable for home-use.

Gene expression profile class prediction using linear Bayesian classifiers.

Due to recent advances in DNA microarray technology, using gene expression profiles, diagnostic category of tissue samples can be predicted with high accuracy. In this study, we discuss shortcomings of some existing gene expression profile classification methods and propose a new approach based on linear Bayesian classifiers. In our approach, we first construct gene-level linear classifiers to identify genes that provide high class-prediction accuracies, i.e., low error rates. After this screening phase, starting with the gene that offers the lowest error rate, we construct a multi-dimensional linear classifier by incorporating next best-performing genes, until the prediction error becomes minimum or 0, if possible. When we compared classification performance of our approach against prediction analysis of microarrays (PAM) and support vector machines (SVM) based approaches, we found that our method outperforms PAM and produces comparable results with SVM. In addition, we observed that the gene selection scheme of PAM could be misleading. Albeit SVM achieves relatively higher prediction performance, it has two major disadvantages: Complexity and lack of insight about important genes. Our intuitive approach offers competing performance and also an efficient means for finding important genes.

Determining a continuous marker for sleep depth.

Detection and quantification of sleep arousals is an important issue, as the frequent arousals are known to reduce the quality of sleep and cause daytime sleepiness. In typical sleep staging, electroencephalograph (EEG) is the core signal and based on the visual inspection of the frequency content of EEG, non-rapid eye movement sleep is staged into four somewhat rough categories. In this study, we aimed at developing a continuous marker based on a more rigorous spectral analysis of EEG to measure or quantify the depth of sleep. In order to develop such a marker, we obtained the time-frequency map of two EEG channels around sleep arousals and identified the frequency bands that show the most change during arousals. We then evaluated classification performance of the potential signals for representing the depth of sleep, using receiver operating characteristic analysis. Our comparisons based on the area under the curve values revealed that the sum of absolute powers in alpha and beta bands is a good continuous marker to represent the depth of sleep. Higher values of this marker indicate low-quality sleep and vice versa. We believe that use of this marker will lead to a better quantification of sleep quality.

Feedback network controls photoreceptor output at the layer of first visual synapses in Drosophila.

At the layer of first visual synapses, information from photoreceptors is processed and transmitted towards the brain. In fly compound eye, output from photoreceptors (R1-R6) that share the same visual field is pooled and transmitted via histaminergic synapses to two classes of interneuron, large monopolar cells (LMCs) and amacrine cells (ACs). The interneurons also feed back to photoreceptor terminals via numerous ligand-gated synapses, yet the significance of these connections has remained a mystery. We investigated the role of feedback synapses by comparing intracellular responses of photoreceptors and LMCs in wild-type Drosophila and in synaptic mutants, to light and current pulses and to naturalistic light stimuli. The recordings were further subjected to rigorous statistical and information-theoretical analysis. We show that the feedback synapses form a negative feedback loop that controls the speed and amplitude of photoreceptor responses and hence the quality of the transmitted signals. These results highlight the benefits of feedback synapses for neural information processing, and suggest that similar coding strategies could be used in other nervous systems.

Segmentation of cDNA microarray spots using markov random field modeling.

MOTIVATION: Spot segmentation is a critical step in microarray gene expression data analysis. Therefore, the performance of segmentation may substantially affect the results of subsequent stages of the analysis, such as the detection of differentially expressed genes. Several methods have been developed to segment microarray spots from the surrounding background. In this study, we have proposed a new approach based on Markov random field (MRF) modeling and tested its performance on simulated and real microarray images against a widely used segmentation method based on Mann-Whitney test adopted by QuantArray software (Boston, MA). Spot addressing was performed using QuantArray. We have also devised a simulation method to generate microarray images with realistic features. Such images can be used as gold standards for the purposes of testing and comparing different segmentation methods, and optimizing segmentation parameters. RESULTS: Experiments on simulated and 14 actual microarray image sets show that the proposed MRF-based segmentation method can detect spot areas and estimate spot intensities with higher accuracy.

Use of meixner functions in estimation of Volterra kernels of nonlinear systems with delay.

Volterra series representation of nonlinear systems is a mathematical analysis tool that has been successfully applied in many areas of biological sciences, especially in the area of modeling of hemodynamic response. In this study, we explored the possibility of using discrete time Meixner basis functions (MBFs) in estimating Volterra kernels of nonlinear systems. The problem of estimation of Volterra kernels can be formulated as a multiple regression problem and solved using least squares estimation. By expanding system kernels with some suitable basis functions, it is possible to reduce the number of parameters to be estimated and obtain better kernel estimates. Thus far, Laguerre basis functions have been widely used in this framework. However, research in signal processing indicates that when the kernels have a slow initial onset or delay, Meixner functions, which can be made to have a slow start, are more suitable in terms of providing a more accurate approximation to the kernels. We, therefore, compared the performance of Meixner functions, in kernel estimation, to that of Laguerre functions in some test cases that we constructed and in a real experimental case where we studied photoreceptor responses of photoreceptor cells of adult fruitflies (Drosophila melanogaster). Our results indicate that when there is a slow initial onset or delay, MBF expansion provides better kernel estimates.

Reliability analysis of microarray data using fuzzy c-means and normal mixture modeling based classification methods.

MOTIVATION: A serious limitation in microarray analysis is the unreliability of the data generated from low signal intensities. Such data may produce erroneous gene expression ratios and cause unnecessary validation or post-analysis follow-up tasks. Therefore, the elimination of unreliable signal intensities will enhance reproducibility and reliability of gene expression ratios produced from microarray data. In this study, we applied fuzzy c-means (FCM) and normal mixture modeling (NMM) based classification methods to separate microarray data into reliable and unreliable signal intensity populations. RESULTS: We compared the results of FCM classification with those of classification based on NMM. Both approaches were validated against reference sets of biological data consisting of only true positives and true negatives. We observed that both methods performed equally well in terms of sensitivity and specificity. Although a comparison of the computation times indicated that the fuzzy approach is computationally more efficient, other considerations support the use of NMM for the reliability analysis of microarray data. AVAILABILITY: The classification approaches described in this paper and sample microarray data are available as Matlab( TM ) (The MathWorks Inc., Natick, MA) programs (mfiles) and text files, respectively, at http://rc.kfshrc.edu.sa/bssc/staff/MusaAsyali/Downloads.asp. The programs can be run/tested on many different computer platforms where Matlab is available. CONTACT: asyali@kfshrc.edu.sa.

Assessment of reliability of microarray data and estimation of signal thresholds using mixture modeling.

DNA microarray is an important tool for the study of gene activities but the resultant data consisting of thousands of points are error-prone. A serious limitation in microarray analysis is the unreliability of the data generated from low signal intensities. Such data may produce erroneous gene expression ratios and cause unnecessary validation or post-analysis follow-up tasks. In this study, we describe an approach based on normal mixture modeling for determining optimal signal intensity thresholds to identify reliable measurements of the microarray elements and subsequently eliminate false expression ratios. We used univariate and bivariate mixture modeling to segregate the microarray data into two classes, low signal intensity and reliable signal intensity populations, and applied Bayesian decision theory to find the optimal signal thresholds. The bivariate analysis approach was found to be more accurate than the univariate approach; both approaches were superior to a conventional method when validated against a reference set of biological data that consisted of true and false gene expression data. Elimination of unreliable signal intensities in microarray data should contribute to the quality of microarray data including reproducibility and reliability of gene expression ratios.

Assessment of closed-loop ventilatory stability in obstructive sleep apnea.

Previous studies on ventilatory control in obstructive sleep apnea (OSA) have generally indicated depressed chemosensitivity, implying greater stability of the chemical control of breathing in these subjects. However, these results were based on tests involving steady-state or quasi-steady measurements obtained in wakefulness. We have developed a method for assessing the dynamic stability characteristics of chemoreflex control in OSA patients during sleep. While continuous positive airway pressure was applied to stabilize the upper airways, acoustically stimulated arousals were used to perturb the respiratory system during sleep. The fluctuations in esophageal pressure that ensued were analyzed, using a closed-loop minimal model, to estimate the chemoreflex loop impulse response (CLIR). Tests using simulated data confirmed the validity of our estimation algorithm. Application of the method to arousal responses measured in six OSA and five normal subjects revealed no statistically significant differences in gain margins and loop gain magnitudes between the two groups. However, the CLIR in the OSA subjects exhibited faster and more oscillatory dynamics. This result implies that, in addition to unstable upper airway mechanics, an underdamped chemoreflex control system may be another important factor that promotes the occurrence of periodic obstructive apneas during sleep.